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SECTION  1 
INTRODUCTION 


The  work  presented  here  continues  the  study  of  the  pilot  Induced 
oscillation  (PIO)  phenomenon  begun  In  Reference  1.  The  results  pre¬ 
sented  In  this  report  for  the  lateral-directional  mode  are  comple¬ 
mentary  to,  but  not  Identical  with,  those  for  longitudinal  dynamics. 

Until  much  more  data  of  quality  can  be  obtained  from  flight  test, 
it  is  probably  not  possible  to  determine  the  validity  of  these 
theories.  In  view  of  the  complexities  associated  with  the  study  of 
lateral-directional  dynamics.  It  would  be  perhaps  premature  to  suggest 
that  this  work  represents  a  mature  theory.  Rather,  It  is  offered  for 
the  assistance  of  others  In  the  field  who  may  have  the  Interest  and 
resources  to  test  it  in  actual  practice. 

The  method  proposed  In  this  report  for  the  assessment  of  lateral- 
directional  PIO  appears  to  apply  equally  well  to  the  longitudinal 
mode.  To  this  author,  this  suggests  that  there  may  be  an  additional 
mode,  not  pursued  In  Reference  1,  by  which  PIO  can  develop — specifi¬ 
cally  the  "single  loop"  tracking  of  pitch  attitude.  It  is  conjectured 
that,  for  this  attitude-only  moda,  normal  acceleration  dynamics  are 
irrelevant,  and  that  this  mode  probably  only  will  be  seen  when  time 
delay  Induces  significant  phase  lag  within  the  bandwidth  of  piloted 
control.  This  delay  might  be  due  to  flight  control  system  nonlineari¬ 
ties,  digital  flight  control  system  delays,  or  equivalent  delay  due  to 
higher  order  system  dynamics. 

The  results  of  this  study  are  very  closely  related  to  methods  for 
PIO  assessment  that  have  existed  for  many  years  (e.g..  Reference  2). 
For  example,  it  has  been  said  (Reference  3)  that  if  "wlde"-band  con¬ 
trol  with  a  pure  gain  pilot  model  is  possible,  then  PIO  is  unlikely. 
The  present  work  supports  that  position  and  provides  a  quantitative 
definition  of  required  system  bandwidth. 
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The  two  principal  contributions  of  this  report  and  Reference  1  to 
the  theory  of  PIO  are; 

1.  To  present  a  theory  that  can  explain  why  the  "primitive" 
pilot  model  form  is  appropriate  for  PIO  study,  and 

2.  To  present  a  quantitative  method  for  the  prediction  of 
pilot-aircraft  systm  bandwidth  sufficient  for  the  assess¬ 
ment  of  PIO  tendencies. 

In  the  world  of  flight  test  applications.  It  is  Important  to 
maintain  an  awareness  of  how  difficult  it  is  to  validate  any  theory 
for  PIO.  There  are  no  test  methods  currently  available  which  can 
guarantee  that  PIO  can  be  found  when  it  is  truly  a  potential  problem, 
or  that  will  ensure  that  PlO-free  configurations  will  be  identified  as 
such.  This  author  is  not  yet  prepared  to  accept  the  notion  currently 
in  vogue  that  any  airplane  can  have  a  PIO  problem.  Any  airplane  can 
be  made  to  oscillate,  certainly. 

It  is  noted  that  there  are  many  other  factors  of  possible  sig¬ 
nificance  to  PIO  in  addition  to  the  aircraft  flight  control  system 
dynamics  addressed  in  this  report.  Most  of  these  are  considered  in 
References  2  and  3.  The  work  of  Reference  1  and  this  report,  however, 
is  directed  toward  a  study  of  aircraft  dynamics  in  order  to  support 
the  development  of  critical  requirements  for  MIL-F-8785C.  These 
generally  are  those  which  could  affect  the  preliminary  design  process 
and  enable  the  development  of  better  tradeoffs  between  aerodynamic 
configuration  and  flight  control  system  requirements. 

The  method  for  the  assessment  of  PIO  tendency  is  presented  for¬ 
mally  in  Section  IV.  The  estimation  of  system  bandwidth  (Section  II) 
is  considered  to  be  crucial  to  the  successful  prediction  of  PIO  (and 
to  handling  qualities,  in  general).  The  discussion  of  the  "primitive" 
pilot  model  given  in  Section  III  is  a  restatement  of  previous  discus¬ 
sions  in  References  1  and  4.  This  material  is  included  here  for  the 
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benefit  of  those  who  may  wish  to  pursue  the  pilot  modeling  problem  In 
future  research  efforts.  Existing  data  are  examined  In  Section  V  to 
evaluate  the  basic  usefulness  of  the  proposed  FIO  assessment  method. 
The  appendix  contains  a  discussion  of  one  of  the  important  handling 
qualities  data  sets.  These  data  are  for  longitudinal  control  In  the 
approach  and  landing  task.  The  discussion  Is  Included  In  this  report 
because  the  data  set  was  derived  from  a  series  of  flight  tests  which 
were  Intended  to  search  for  PIO.  These  data  were  important  to 
refining  the  estimation  of  closed  loop  bandwidth. 
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SECTION  11 

ESTIMATION  OF  PILOT-AIRCRAFT  SYSTEM  BANDWIDTH 


The  bandwidth  of  the  pilot-aircraft  system  Is  the  single  most 
Important  parameter  to  th~i  success  of  handling  qualities  prediction 
methods  and  to  the  quantitative  analysis  of  PIO.  The  PIO  theory  of 
Reference  1  offered  no  convenient  method  for  bandwidth  prediction.  It 
was  suggested  that.  In  order  to  apply  the  theory,  available  methods 
for  pilot-vehicle  system  analysis  could  be  used  to  estimate  whether 
closed  loop  resonance  was  likely,  and  the  resonant  frequency.  Because 
of  the  difficulty  of  accomplishing  this  In  practice,  the  development 
of  a  recommended  PIO  design  specification  (published  as  an  appendix  In 
Reference  5)  left  the  definition  of  system  bandwidth  purposely  vague. 

Reference  5  did  offer  a  method  for  the  estimation  of  pilot- 
vehicle  system  bandwidth  In  non-PIO  conditions.  This  was  the  key  to 
the  method  proposed  for  handling  qualities  assessment  and  prediction 
In  Category  A  tracking  tasks  with  Type  IV  aircraft.  The  method  pro¬ 
posed  was  based  on  tiro  observations: 

1.  During  the  process  of  selecting  appropriate  pilot  models  for 
the  analysis  of  pre-PIO  tracking  dynamics  In  Reference  1, 
the  Bode  amplitude  properties  seemed  to  dominate  the  selec¬ 
tion  process;  aircraft  phase  response  appeared  to  be  a 
secondary  consideration.  (It  Is  recalled  that  this  process 
Is  nothing  more  than  the  application  of  the  "adjustment 
rules"  from  Reference  6  to  estimate  realistic  parameters  for 
the  servo  model  description  of  the  human  pilot.) 

2.  The  crossover  frequency  (bandwidth)  data  summarized  In 
Reference  6  for  single  axis  tracking  suggest  that  crossover 
frequency  Is  almost  completely  parameterized  by  the  Bode 
amplitude  properties  of  the  controlled  element  (aircraft) 
dynamics.  Whether  this  observation  had  been  made  by  others 
Is  unknown;  It  was  a  revelation  to  this  author. 
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Work  subsequent  to  that  of  Reference  5  suggested  that  the  cross¬ 
over  frequency  predicted  by  the  formula  of  Reference  5  was  a  good 
first  approximation  to  the  PIO  frequency — particularly  when  PIO  devel¬ 
oped  as  the  result  of  linear  system  dynamics.  While  nonllnearltles  In 
the  flight  control  system,  such  as  saturation,  might  affect  the  fully 
developed  PIO  frequency,  the  possibility  exists  that  the  Initial  PIO 
frequency  may  be  equal  to  the  crossover  frequncy.  In  the  development 
of  design  specifications,  the  properties  of  fully  developed  PIO  (e.g., 
frequency  and  amplitude  of  the  motion)  are  of  no  real  Interest.  Only 
those  conditions  which  precipitate  PIO  are  of  value  to  the  development 
of  criteria  which  can  eliminate  PIO  as  a  problem.  Thus,  the  crossover 
frequency,  computed  by  the  method  of  Reference  5,  can  be  used  as  the 
resonant  frequency  for  evaluation  of  the  proposed  PIO  design 
criterion. 

A.  THE  CROSSOVER  FREQUENCY  FORMULA 

The  basic  relation  for  crossover  frequency,  as  a  function  of 
aircraft  dynamics  was  derived  In  Reference  5  based  on  measurements  of 
crossover  frequency  published  In  Reference  6.  For  pitch  attitude 
control.  It  was  determined  that  Is  parameterized  by  the  average 
slope  of  the  pitch  attitude  Bode  amplitude  response  over  the  frequency 
Interval  of  2  to  6  rad Ians/ second.  A  limit,  equal  to  6  radians/second 
was  Imposed  on  u^. 

Evaluation  of  the  crossover  frequency  formula  against  other  data 
sources  suggests  that  the  formula  should  be  slightly  modified  as 
follows: 

1.  The  slope  of  |0/Fgg(ju)|  should  be  averaged  over  1  to 
6  rad Ians/ second. 

2.  The  restriction  on  the  maximum  value  of  u)  should  be 

c 

removed. 


■"•'T 


5 


The  first  modification  brings  the  formula's  predictions  more  in  line 
with  observations  based  on  the  LAHOS  data  (discussed  in  the  appendix) . 
The  second  modification  was  suggested  based  on  discussions  with  Air 
Force  Flight  Test  personnel  regarding  F-15  flight  test  data.  The 
resulting  formula  Is: 

,  rad Ians/ second 

where  S  -  average  slope  of  |0/Fgg(jft)) | ,  In  decibels  per  octave,  on  the 
region  1  <  u  <  6  rad Ians/ second. 

The  data  of  Reference  6  on  which  the  formula  Is  based  were  from  a 
lateral  control  task,  not  unlike  a  roll  control  task.  It  Is  expected, 
therefore,  that  the  formula  will  be  suitable  for  applications  to 
lateral-directional  control  problems.  For  lateral-directional  con¬ 
trol,  the  slope,  S,  will  be  the  average  slope  of  | i|'/F^(ju)  | . 

Throughout  the  remainder  of  this  report,  the  slope,  S, 
is  computed  as  follows: 

S  =  j  [a(2)  -  A(l)  +  A(3)  -  A(1.5)  +  A(4)  -  A(2) 

+  A(5)  -  A(2.5)  +  A(6)  -  A(3)] 

where 

A»  in  decibels,  or  |4'/F^(ju))  | . 

B.  CONFIRMATION  OF  THE  CROSSOVER  FREQUENCY  FORMULA 

There  are  very  little  data  available  which  are  suitable  for 
establishing  the  validity  of  the  above  formula  for  the  roll  tracking 
task.  Two  restricted  data  sets  will  be  considered. 


In  Reference  7,  Durand  and  Jex  published  estimates  of  the  domi¬ 
nant  closed  loop  natural  frequency  observed  over  short  time  Intervals 
In  a  roll  tracking  experiment.  Their  values  should  approximate  the 
crossover  frequency  and  can  be  used  to  test  the  formula.  The  simu¬ 
lated  roll  dynamics  were: 


*^AS 


(s) 


s  (s  +  ~) 


Three  values  of  the  roll  subsidence  time  constant  were  tested. 
For  one  of  these  the  pilot  subject  A  was  asked  to  track  "normally," 
then  Intentionally  use  a  hlgh-galn  technique,  with  a  final  reversion 
to  the  normal  tracking  technique.  The  average  closed  loop  frequency 
observed  In  each  case,  Is  compared  In  the  table  below  with  the 

crossover  frequency  predicted  using  the  above  formula. 


Pilot  Control  Technique  (Pilot  A) 


Revert-to- 

Normal 

Normal 

Hi-Gain 

Normal 

Normal 

1 

^R 

1.0 

2.5 

2.5 

2.5 

10.0 

"cL 

2.5 

3.0 

3.9 

3.1 

4.1 

U) 

c 

3.34 

3.8 

3.8 

3.8 

4.43 

It  appears  that  by  formula.  Is  generally  larger  than  actual  closed 
loop  oscillation  frequency,  except  for  the  worst  case,  hlgh-galn 
condition  where  the  agreement  Is  very  good. 

In  general,  probably  depends  on  the  pilot.  Durand  and  Jex 
showed  that  systematic  variations  In  occur,  depending  upon  the 


7 


pilot.  The  table  below  shows  obtained  for  each  of  four  pilots. 
The  average  pilot  rating,  R,  Is  also  shorn  In  each  case  (modified 
Cooper  scale) : 


Pilot 

R 

“CL 

(rad/ sec) 

Predicted 

R 

Range 

Median 

D 

2.25 

2. 3-3. 4 

2.8 

3.65 

A 

2.35 

3. 0-3. 2 

3.1 

3.75 

B 

2.50 

3. 1-3.8 

3.4 

4.15 

C 

3.75 

4. 0-4.4 

4.3 

4.75 

The  overall  average  ■  3.4  Is  In  very  good  agreement  with  the  id 
formula.  In  the  table,  the  data  are  ordered  according  to  pilot 
rating.  Mote  that  rating  Increases  monotonlcally  with  frequency  of 
loop  closure.  This  supports  the  theory  of  Reference  5  for  rating 
prediction  (which  concludes  that  rating  Increases  with  phase  lag  of 
Che  controlled  element,  measured  at  the  crossover  frequency).  This 
variation  Is  shown  In  Figure  1. 

The  ratings  that  would  be  predicted  for  each  of  these  cases, 
using  u  •  and  the  pilot  rating  function  from  Reference  5  for 
pitch  tracking,  are  shown  In  the  table.  These  predicted  ratings  do 
not  agree  with  the  actual  ratings.  Based  upon  limited  data.  It 
appears  that  Che  ratings  received  In  roll  tracking  are  not  so  sensi¬ 
tive  to  controlled  element  phase  as  In  Che  case  of  pitch  control.  Use 
of  the  rating  function  from  Reference  5  for  the  prediction  of  roll 
handling  qualities  seems  to  yield  conservative  predictions.  However, 
the  Durand-Jex  data  do  support  Che  crossover  frequency  formula  for 
applications  to  roll  tracking. 

An  Interesting,  brief  flight  test  of  lateral-directional  handling 
qualities  was  conducted  as  part  of  a  larger  test  and  reported  In 
Reference  8.  Smith  et  al.  analyzed  these  data  and  supplied  the  author 
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DOMINANT  CLOSED  LOOP  FREQUENCY 
WqL  RAD/SEC 


Figure  1.  Closed  Loop  Frequency  vs.  Pilot  Rating 

with  some  of  their  preliminary  results.  These  Included  response  time 
histories  and  power  spectral  density  (PSD)  plots  of  lateral  stick 
force  for  four  configurations  with  which  PIO  was  experienced  in 
landing.  These  were  roll-only  PlOs  to  a  good  approximation;  the  dutch 
roll  mode  was  canceled  In  the  transfer  function. 

The  following  table  summarizes  the  pilot  ratings,  the  PIO  fre¬ 
quency  (from  the  PSDs  or  time  traces),  and  the  crossover  frequency 
determined  by  application  of  the  formula  given  in  Section  IIA: 
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1 

j 


Conf ig . 

Cooper- 

Harper 

Rating 

PIO  Frequency 

(rad/sec) 

Flight  Test 

Formula 

L8 

5/6/9 

2.3 

2.38 

LllC 

6/9 

3. 5-4.0* 

3.58 

LlAB 

5/7/10 

1.7 

1.99 

L16A 

8 

3.5 

3.22 

*Estimated  for  pilot  A,  record  15  where  PIO  was  a  definite 
problem. 


It  appears  that  these  data  support  the  formula  proposed  for  the  esti¬ 
mation  of  system  bandwidth. 

One  other  case  is  worth  noting  since  it  isn't  generally  acces¬ 
sible  in  the  open  literature.  For  the  Space  Shuttle  Orbiter,  a  mild 
longitudinal  PIO  was  experienced  on  the  fifth  free-fllght  in  the 
post-flare  region  of  landing.  Based  on  power  spectral  density  data 
for  free-f light  5,  the  PIO  frequency  (pitch  mode)  was  about  3.5  radians/ 
second.  Using  the  Orbiter  dynamics  and  the  crossover  frequency 
formula,  it  can  be  shown  that  ^  3.8.  The  agreement  with  flight 
test  is,  therefore,  very  good. 

It  is  noted  that  the  bandwidth  formula  indicates  that  pilot- 
aircraft  system  bandwidth  will  not  be  less  in  approach  and  landing 
than  in,  for  example.  Category  A  flight  phases.  This  is  an  Important 
result  which  should  have  significant  Impact  on  the  character  of  design 
specifications  for  handling  qualities  in  the  landing  task. 

Other  aircraft  examples  will  be  shown  in  Section  V  which  further 
support  the  crossover  frequency  formula. 
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C.  THEORETICAL  MODELS  FOR  CROSSOVER  FREQUENCY 

The  formula  above  for  crossover  frequency  appears  to  be  generally 
adequate  for  studies  of  aircraft  handling  qualities  In  attitude  con¬ 
trol.  Despite  Its  success  with  the  limited  testing  It  has  received  to 
date,  the  formula  remains  empirical. 

To  have  genuine,  widespread  value  to  handling  qualities  analysis. 
It  Is  Important  that  the  technique  used  for  bandwidth  prediction  be 
applied  across  all  conditions  of  task,  aircraft  dynamics,  and  external 
Influences.  This  will  require  the  development  of  a  generalized  model 
for  closed  loop  bandwidth.  This  goal  may  be  beyond  our  reach  at  the 
present  time.  Still,  for  purposes  of  discussion  It  may  be  of  value  to 
consider  a  very  limited  aspect  of  the  theoretical  problem. 

In  Reference  4  It  was  hypothesized  that  a  pilot-aircraft  model  of 
the  form  shown  In  Figure  2  could  possibly  account  for  most  of  the 
trends  found  In  the  human  pilot  data  base  provided  by  Reference  6.  It 
was  speculated  that  If  appropriate  nonlinear It les  were  Included  In  the 
model— particularly  In  the  rate  control  path — the  model  would  apply, 
without  parameter  variations,  over  the  entire  range  of  controlled 
elements  tested  In  Reference  6.  These,  however,  were  never  found 
(although  some  tantalizing  results  were  produced). 

The  version  of  the  model  shown  In  Figure  2  Is  Intended  to  repre¬ 
sent  a  linearization  of  the  more  complete,  but  abstract,  model.  When 
the  Inner  (rate)  loop  Is  closed,  the  block  diagram  structure  Is  Iden¬ 
tical  to  the  single  loop,  unity  gain  feedback  model  generally  assumed 
for  the  servo  model  description  of  the  pilot-vehicle  system.  The 
resulting  open  loop  transfer  function  [r/r^(s)]  can  be  called 

Y  (s)  Y  (s),  as  Is  done  with  the  servo  description.  When  this  Is 
P  ^ 

done,  then  ^^(s)  becomes  a  function  of  Y^(s),  the  controlled  element 
dynamics.  The  claim  was  made  In  Reference  4  that  this  functional 
behavior  can  account  for  the  requirements  Imposed  on  parameterization 
of  the  servo  model  by  the  "adjustment  rules"  of  Reference  6. 
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Y  (s)  =  controlled  element  transfer  function 
c 


-.095s 


G_(s)  = 


VJ  VO/  —  y  i  —  y  ^  »  1&CU1.V/U1U3 

"  (fo  - ')  [(Sf  - 

=  equivalent  rate  washout  time  constant 


>  neuromuscular  system 


Ti  =  time  delay,  rate  path 
T2  =  time  delay,  position  path 

^1>  ^2  **  pilot  gain  constants,  rate  and  position  paths, 
respectively 


Figure  2.  First  Order,  Multiple  Loop  Model  for 
Pilot-Aircraft  System  Dynamics 
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The  purpose  of  discussing  the  multiple  loop  pilot  model  theory  In 
this  report  Is  to  demonstrate  that  there  Is  probably  a  systematic, 
analytical  technique  that  will  enable  the  accurate  a  priori  estimation 
of  system  bandwidth,  ui^,  for  conditions  of  general  handling  qualities 
Interest. 

Assume  the  following  values  for  the  pilot  model  parameters  of 
Figure  2; 


Ki  =«  1.0  K2  -  7.0  1/T^  -  .05 

Ti  =  .1  T2  ■  .2 

With  these  values  fixed,  the  model  of  Figure  2  can  be  used  to  compute 
the  open  loop  system's  frequency  response  for  any  controlled  element. 
Figures  3a,  3b,  and  3c  show  the  results  for  Y^(s)  ■  K^,  K^/s,  and 
K^/s^.  Representative  measured  values  of  YpCjiu)  Y^(jw)  from  Refer¬ 
ence  6  are  shown  for  comparison.  For  clarity,  the  data  variances  from 
Reference  6  are  not  shown. 

It  is  apparent  that  the  prediction  of  crossover  frequency  for  all 
three  configurations  is  very  good.  The  correlation  between  model  and 
experiment  is  generally  very  good  for  Y^Cs)  *  K^,  and  for  w  <  o)^, 

K^/s.  For  Y^(s)  *  K^/s^,  the  phase  angle  correlation  Is  poor.  Better 
agreement  could  have  been  obtained  by  modifying  the  neuromuscular 
system  parameters. 

Unfortunately,  the  same  model  parameters  do  not  work  as  well  for 

the  other  Y  (s)  tested  in  Reference  6  (tabulated  on  the  following 
c 

page).  This  Is  to  be  expected,  since  the  model  Is  but  a  crude  approx¬ 
imation  to  that  required  by  hypothesis. 

Nevertheless,  It  Is  Instructive  to  visualize  how  such  a  model 
might  be  used  to  explore  handing  qualities  sensitivity  to  pilot  tech¬ 
nique.  In  the  following  table,  the  change  In  crossover  frequency  due 
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to  a  50  percent  decrease  In  the  pilot's  rate  gain  is  shown  as  a  func¬ 
tion  of  Y^(s).  All  other  pilot  model  parameters  are  constant  and 
equal  to  those  used  In  Figures  3a,  3b,  and  3c. 


50  Percent  Decrease  In 


Y^(s) 

Percent  Change  in 

K 

2332 

increase 

c 

K  /s 

102 

increase 

c 

K  /s2 

3% 

decrease 

K  /(s-2) 
c 

402 

increase 

K  /s(s-1.5) 
c 

13% 

decrease 

K^/(s-.425) 

02 

change 

If  it  were  known  how  pilot  parameters  such  as  varied  as  a 
function  of  stress,  g-loadlng,  etc.,  then  a  general  model  for  the 
estimation  of  would  permit  prediction  of  the  resulting  effect  on 
handling  qualities.  This  assumes  that  a  parameter  such  as 
y  ♦/F^gCjUc)  has  a  one-to-one  correspondence  with  pilot  rating. 

It  Is  not  sufficient  to  merely  assume  a  specific  change  In  and 
see  what  the  effect  on  handling  qualities  would  be.  There  are  at 
least  two  sensitivities  that  must  be  jointly  considered: 

•  the  sensitivity  of  pilot  model  parameters  to  factors  such  as 
stress,  muscle  tension,  external  constraints,  etc.  (which 
invoke  physiological  considerations) ,  and 

•  the  sensitivity  of  system  bandwidth  to  variations  In  pilot 
parameters  (which  Is  a  function  oZ  aircraft  dynamics). 
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Figure  3a.  Open  Loop  Pilot-Vehicle  System  Frequency 

Response:  Y  (s)  •  K 
c  c 


^  Yp(i<i»)Yc(i«)~DEGS  |Yp{jai)Yc(j«)|~DB 


T  . 


RAD/SEC 


Figure  3b.  Open  Loop  Pilot-Vehicle  System  Frequency 
Response;  Y^(s)  ■  K^/s 
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tYp(i«)Yc{iw)~DE6  |Yp(i«)Yc(j«)| 


As  a  final  conment,  the  sensitivity  data  tabulated  previously 
support  the  qualitative  arguments  made  in  Reference  5  for  the  first 
three  controlled  elements.  That  Is,  o)^  Is  extremely  sensitive  to  rate 
gain  when  Y^(s)  -  K^,  not  very  sensitive  for  K^/s,  and  almost  Indepen¬ 
dent  of  rate  gain  for  K^/s^. 


SECTION  III 
THE  "PRIMITIVE  PILOT" 


The  pure  gain  pilot  model,  often  suggested  and  used  In  the  PIO 
literature,  was  originally  proposed  because  It  seemed  to  be  consistent 
with  limited  PIO  case  history  data.  Its  origins.  In  other  words,  are 
empirical.  Within  the  general  context  of  the  servo  theory  for  human 
dynamics,  the  pure  gain  model  appears  to  be  out  of  place.  There  Is, 
for  example,  no  obvious  way  that  It  can  be  derived  from  the  genral 
servo  model  through  application  of  the  adjustment  rules  of  Refer¬ 
ence  6.  This  connection  between  theory  and  practical  observation  Is 
an  Important  one  to  make.  It  Is  reasonable  to  expect  that  the 
degeneration  of  pilot  model  form  should  somehow  reflect  the  physical 
mechanisms  by  which  PIO  can  be  Initiated.  If  these  can  be  understood, 
then  there  Is  hope  that  PIO  can  be  eliminated  from  future  designs. 

One  theoretical  explanation  which  justifies  the  pure  gain  model 
for  fully  developed  PIO  was  offered  In  Reference  1.  It  was  noted 
there  that  the  multiple  loop  pilot  model  (outlined  In  the  previous 
section)  appears  to  be  structurally  consistent  with  plausible  PIO 
physical  behavior  and  with  the  accompanying  degeneration  of  pilot 
model  form.  In  the  multiple  loop  model,  each  mode  of  feedback  control 
derives  from  an  Independent  physiological  sensor  element.  Thus, 
visually  derived  error  rate  Is  not  treated  as  the  time  derivative  of 
error  since  there  exist  different  sensor  dynamics,  delays,  and  cere¬ 
bral  pathways.  In  the  Incipient  PIO  state  If  the  pilot  simply 
switches  off  all  feedbacks  except  that  which,  to  him,  is  most  task 
critical,  then  the  resulting  pilot  model  Is 

Y  (s)  -  K  e"^® 

P  P 

In  other  words,  the  "primitive  pilot"  model  may  be  merely  the  forward- 
path  pilot  transfer  function  of  the  one  remaining  feedback.  The 
concept  of  the  primitive  control  Is  a  cornerstone  of  the  theory  for 
longitudinal  PIO  In  Reference  1  (where  the  remaining  single  loop  was 

I 
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postulated  to  be  normal  acceleration).  In  Reference  1,  available  data 
suggested  a  nonzero  time  delay  should  be  retained.  It  Isn't  clear 
whether  this  will  also  be  true  when  the  FIO  Is  due  to  single  axis 
pitch  or  roll  tracking.  Throughout  the  remainder  of  this  report  x  *  0 
will  be  assumed,  for  viant  of  better  information. 

Following  the  postulates  of  Reference  1,  It  Is  assumed  that  the 
pilot  regression  to  the  primitive  form  of  single  mode  control  follows 
the  development  of  substantial  aircraft  resonance  In  attitude 
response.  This  can  result  from  open  or  closed  loop  Inputs. 


SECTION  IV 

A  PROPOSED  METHOD  FOR  PIO  ASSESSMENT 


As  noted  in  the  Introduction,  the  method  proposed  here  for  PIO 
assessment  Is  not  startlingly  different  from  traditional  methods 
(e.g..  Reference  2).  The  essential  difference  Is  that  a  quantitative 
definition  of  required  system  bandwidth  Is  offered.  Also,  the  pro¬ 
posed  assessment  method  attempts  to  establish  a  positive  relation 
between  the  essential  behavior  of  the  real  (nonlinear)  flight  control 
system  and  the  Idealized  (linear)  one. 

The  method  proposed  la  confirmed  In  Section  V  against  available 
case  history  data  for  lateral-directional  PIO.  There  Is  no  obvious 
reason  to  preclude  the  application  to  longitudinal  dynamics,  as  well. 

In  the  case  of  longitudinal  PIO,  it  Is  suggested,  but  not  proven, 
that  PIO  In  pitch-only  tracking  Is  possible  according  to  the  proposed 
theory.  It  appears  that  this  would  require  sufficient  time  delay 
(equivalent,  due  to  higher  order  dynamics,  or  real)  to  make  the 

0  single  loop  unstable  at  m  (from  the  formula) .  It  Is  hypothe- 

Ciw  c 

sized  that  this  would  be  true  regardless  of  normal  acceleration  dynam¬ 
ics;  If  so.  It  would  constitute  another  PIO  mode  distinct  from  that 
given  by  the  theory  of  Reference  1.  If  this  Is  true,  then  It  should 
be  possible  to  excite  longitudinal  PIO  In  fixed-base  simulation  pro¬ 
vided  sufficient  delay  exists  in  the  0/F_„  transfer  function.  This 

£>0 

contradicts  the  author's  conjecture  In  Reference  1  regarding  PIO  in 
fixed-base  simulation.  In  all  other  respects.  It  Is  believed  that  the 
present  theory  complements  that  of  Reference  1. 

It  should  be  noted  that.  In  the  original  approach  to  the  problem 
of  understanding  and  predicting  lateral-directional  PIO,  the  author 
attempted  to  apply  the  longitudinal  theory  (Reference  1)  with  ())  and 
ay^  substituting  for  6  and  anp>  respectively.  For  the  data  analyzed 
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in  this  study,  It  appears  that  this  approach  Is  not  satisfactory. 

Even  for  the  M2-F2  PIO — which  Involves  substantial  dutch  roll — the 
dynamics  appeared  to  add  nothing  to  understanding  the  problem.  In  all 
cases,  the  (p  ■*  F  loop  dynamics  were  sufficient  for  PIO  evaluation. 
The  use  of  rudder  pedal  contributed  to  control  difficulties  with  the 
M2-F2;  It  does  not  appear  to  have  been  the  reason  for  the  airplane's 
PIO  encounters.  Based  upon  undocumented  PIO  experiences  with  the  B-52 
and  KC-135,  this  may  be  the  usual  Involvement  of  rudder  control  In  PIO 
when  the  controls  are  of  the  classical  sort.  We  should,  however,  be 
prepared  for  "rudder"-only  PIO  when  direct  force  control  systems  are 
Introduced. 

A.  BACKGROUND 

A  Cooper-Harper  scale  rating  of  8  represents  a  boundary  between 

(a)  controllability  not  In  question  (R  -  7) ,  and 

(b)  controllability  la  conditionally  dependent  on  pilot  compen¬ 
sation  (R  -  9) . 

It  Is  reasonable  to  expect  that  a  rating  of  8  might  represent  a 
threshold  to  PIO  in  those  cases  where  the  problem  exists. 

The  relation  derived  In  Reference  5  between  Cooper-Harper  rating 
and  ):  6/F  (jo)  )  Is  shown  In  Figure  4.  Note  the  knee  In  the  relation 
at  R  >  8.  The  corresponding  6-phase  Is  very  close  to  -180  degrees. 

This  suggests  that  a  simple  necessary  condition  for  attitude-only 
PIO,  pitch,  or  roll  would  be  that: 

)  or  y  (J“  )  <  -180  degrees 
^ES  ^AS  "" 
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Figure  4.  Cooper-Harper  Rating  vs.  ):  0/F_„  (jiu  ) 
(from  Reference  5)  ^ 


for  linear  system  dynamics.  Thus,  one  can  use  a  pure  gain  pilot  model 
and  look  for  locus  crossings  of  the  imaginary  axis  on  a  root  locus 
plot.  When  the  frequncy  of  axis  crossing  is  less  than  — by  formula — 

then  the  crossing  conditions  represent  a  potential  PIO  state.  In 
order  to  get  PIO,  something  has  to  excite  the  pilot  to  adopt  the 
primitive  model  form.  The  sketch  following  Illustrates  the  relation 
between  the  frequency  of  axis  crossing,  and  PIO  In  a  roll  tracking 
task  without  dutch  roll.  The  effect  of  increasing  time  delay  Is 
illustrated.  This  could  originate  from  within  the  flight  control 
system  and  be  an  actual  or  equivalent  delay. 


Noc  that  is  a  function  only  of  l8/Fgg(J«))  or  ) 4i/F^ (joi)  |  which 
Is  unaffected  by  time  delay.  Thus,  ui  is  a  constant  for  all  loci  shown. 


CONST 


B.  THE  METHOD 

It  Is  possible  to  generalize  the  above  remarks  to  Include  non¬ 
linear  system  dynamics.  The  resulting  assessment  method  Is  best 
described  by  the  flow  chart  of  Figure  5. 

Previously,  the  discussions  of  system  dynamics  have  been 
restricted  to  closed  loop  dynamics.  It  Is  suggested  that  resonance  In 
roll  (or  pitch)  attitude  due  to  open  loop  control  or  disturbance 
Inputs,  or  to  the  use  of  secondary  controls  (e.g.,  rudder),  must  also 
be  considered  In  a  general  PIO  analysis,  since  this  can  lead  to  PIO 
Initiation.  For  resonance  Induced  by  other  than  the  $  F^  closed 
loop.  It  Is  proposed  that 

1.  When  the  damping  ratio  of  the  dominant  mode,  >  0.2,  then 
there  Is  Insufficient  resonance  to  produce  PIO. 
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ESTIMATE  CRITERION  FREQUENCY, 
FOR  PK)  ASSESSMENT 

•  CLOSED  LOOP  DYNAMICS 

•  OPBI  LOOP  DYNAMICS  S  .2?) 

fYES;  Wc-OiR 
I  NO;  IGNORE 

•  SYSTEM  RESONANCE  INDUCED  BY 
SECONDARY  CONTROLS  (^p  i  .27} 


PERFORM  NONUNEAR  SYSTEMS  ANALYSIS: 

•  UMIT  CYCLE  FREQUENCY  &  AMmjTUDE 

•  PLAUSIBIUTY  OF  LIMIT  CYCLE  ENTRY 


PRIMITIVE"  PILOT  MODEL: 
Yp  (8)  =  Kp 


PK)  POSSIBLE  WITH  INITIAL 
FREQUB^ICY  =  Wg  4  FUaY  - 
DEVELOPED  FREQUENCY  4 
AMPLITUDE  DETmMINED  BY 
LIMIT  CYCLE  BEHAVIOR 
{NOTE  1) 


Figure  5,  PIO  Analysis  Flow  Chart 
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2. 


When  ^  0.2,  then  the  criterion  frequency  for  PIO  analysis 
Is  equal  to  the  resonant  mode  natural  frequency:  '■ 

The  limit  cycle  analysis  necessary  to  determine  limit  cycle 
frequencies,  ,  and  amplitudes  may  be  accomplished  In  any  manner 

Ij 

desired.  With  very  complex  flight  control  systems,  this  will  probably 
require  a  fairly  sophisticated  computer  model  for  the  aerodynamics 
plus  control  system. 

There  Is  necessarily  room  for  art  In  applications  of  the  method 
charted  In  Figure  5.  It  is,  for  example,  all  but  impossible  to  be 
certain  that  limit  cycles  cannot  occur  for  some  Input.  Still,  experi¬ 
ence  suggests  that  past  PIO  occurrences  In  which  control  saturation 
was  a  problem  all  could  have  been  avoided  If  procedures  such  as  that 
of  Figure  5  had  been  used  at  the  time  of  system  design. 

The  following  two  notes  apply  to  the  figure: 

Note  It 

The  potentially  catastrophic  PIO  which  precipitates  structural 
failure  or  pilot  Incapacitation  would  probably  occur  for  this 
case  when  the  system  limits  are  exceeded  prior  to  limiting  due  to 
flight  control  system  (FCS)  saturation  or  aerodynamic  limits. 

The  elimination  of  these  PIOs  will  require  fundamental  revisions 
to  the  FCS  or  to  the  aerodynamic  design. 


Note  2: 


Typically,  expect  that  PIO  resulting  from  this  path  will  be  large 
amplitude,  low  frequency  In  nature,  and  produce  serious  degrada¬ 
tion  of  flight  path  stability.  The  elimination  of  these  PIOs 
will  generally  require  detailed  revision  to  the  FCS;  In  particu¬ 
lar,  saturation  effects  (rate  or  position)  must  be  modified. 


26 


c. 


A  CONJECTURE 


It  is  probably  Important  that  the  possibility  of  PIO  due  to  use 
of  direct  sldeforce  control  (DSFC  or  DFC)  be  considered  within  the 
context  of  Figure  5.  This  author  knows  of  no  PIOs  of  this  nature  that 
have  been  experienced  in  flight  or  simulation. 

As  an  interim  standard  to  prevent  the  occurrence  of  such  PIO 
modes,  it  is  suggested  that  a  modification  of  the  longitudinal  PIO 
theory  be  used  as  follows: 

1.  Determine  the  principal  outer  loop  cue  for  DSFC  use.  This 
may  be  the  offset  between  a  probe  and  drogue  during  aerial 
refueling,  a  HUD-derlved  cue  for  a  ground  attack  mode,  etc. 
Call  this  cue  x. 

2.  Derive  w  based  on  x/F,^  _(s)  where  F  is  the  pilot's  force 

C  Urw  DrL 

input  to  the  direct  force  controller. 

3.  When  ):  ^  -180  degrees,  consider  that  direct 

sldeforce  control  probably  won't  produce  PIO  (barring 
serious  nonlinearities) . 

4.  When  ^  */*’dfC^^‘*’c^  ^  degrees,  then  if  y  ayp/Fpp^Cju)^) 

<  -180  degrees,  conclude  that  PIO  is  a  possibility;  other¬ 
wise,  there  is  probably  not  a  PIO  problem. 
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A. 


SECTION  V 

CONFIRMATION  WITH  AVAILABLE  CASE  HISTORY  DATA 


In  this  section  the  PIO  assessment  methodology  offered  In 
Section  IV  will  be  tested  for  specific  aircraft,  aircraft  dynamics 
simulated  In  flight,  and  for  one  fixed-base  simulation.  Only  the 
lateral-directional  modes  are  considered. 

Complete  nonlinear  systems  analyses  will  not  be  accomplished  as 
part  of  this  validation  study.  Detailed  models  would  be  required  for 
critical  system  nonlinear It les;  generally,  the  data  required  to  con¬ 
struct  these  and  confirm  the  resulting  models  are  not  available  In  the 
literature.  In  those  cases  where  flight  control  system  nonlinear It les 
play  a  significant  role  in  the  development  of  PIO,  other  literature  is 
consulted  or  a  simplified  treatment  Is  used  to  Illustrate  the  PIO 
assessment  methodology. 

A.  YF-16  FIRST  FLIGHT 

Reference  9  contains  an  analysis  of  the  first  flight  PIO  of  the 
YF-16.  The  YF-16  PIO  was  dramatic,  sudden,  and  violent.  It  consisted 
primarily  of  large  rolling  motions  with  fully  saturated  controls.  The 
PIO  occurred  during  what  was  supposed  to  have  been  a  high  speed  taxi 
test  to  evaluate  lateral  control  effectiveness.  An  overspeed  condi¬ 
tion  occurred  and  the  aircraft  became  airborne,  followed  shortly  by 
uncontrolled  PIO.  As’ a  safety  measure  the  pilot  elected  to  add  power, 
gain  altitude,  and  fly,  rather  than  attempt  Immediate  touchdown. 

During  the  PIO,  light  ground  contact  was  made  between  the  right  hori¬ 
zontal  tall  tip  and  the  left  dummy  AIM-9.  The  PIO  time  history  is 
given  In  Figure  6.  In  addition  to  the  data  published  In  Reference  9, 
John  Smith  supplied  to  the  author  the  ^/F.„  frequency  response 
(take-off  condition)  shown  in  Figure  7. 

From  Figure  6  It  appears  that  the  frequency  of  wing-rock  In  the 
preliftoff  period  Is  nearly  5.0  rad Ians/ second.  This  was  conceivably  a 
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factor  in  PIO  Initiation.  It  is  presumed  that  the  wing-rock  resulted 
from  open  loop  control  inputs  from  the  pilot. 

At  main  gear  lift-off  followed  by  left  wing  drop,  aileron  stick 
force  F^g  is  approximately  proportional  to  iJi(t-T).  Following  control 
surface  rate  and  position  saturation,  F^  becomes  approximately  pro¬ 
portional  to  roll  rate  with  no,  or  small,  delay.  This  shift  in  appar¬ 
ent  feedback  cue  architecture  occurs  at  about  t  ■  4  seconds,  aircraft 
tape  time.  For  4  <  t  <  19,  remains  approximately  proportional 

to  p(t) . 

In  the  linear  region  of  the  flight  control  system  and  following 
relaxation  of  stick  force  in  the  post-PIO  region,  the  dominant  closed 
loop  frequency  is  nearly  1  Hz.  In  the  fully  developed  PIO,  there  are 
two  observable  frequencies: 

1.  When  the  control  surfaces  do  not  appear  to  be  position 
saturated,  the  PIO  frequency  is  approximately  5  radians/ 
second. 

2.  When  the  control  surfaces  are  position  saturated,  the  fre¬ 
quency  is  approximately  3.6  rad Ians/ second. 

Sudden  jumps  in  closed  loop  frequency  of  the  types  described  are 
probably  typical  of  large  amplitude,  nonlinear  motions  with  full 
authority  fly-by-wlre  control  systems.  However,  this  will  also  occur 
with  conventional  stability  augmentation  systems  (SAS)  when  limits  of 
rate  or  position  are  exceeded;  the  YF-12  PIO  is  a  good  example  of  this 
(Reference  1) . 

The  crossover  frequency  may  be  determined  by  formula  from  the 
(|)/F^  Bode  plot  of  Figure  7; 
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(ii  *4.18 
c 

y-  (ju  )  *  -186  degrees 

Then  following  the  PIO  assessment  flow  chart  (Figure  5) ,  It  Is  seen 
that  the  (p  closed  loop  Is  potentially  unstable.  The  fully 

developed  PIO  frequency  and  amplitude  are,  therefore,  determined  by 
limit  cycle  states.  Reference  9  Indicates  that  stable  limit  cycle 
frequencies  3.3  £  £  4  are  possible  with  the  pure  gain  pilot  model. 

Since  y  p/F  (ju)  )  is  close  to  the  boundary  value  of  -180  degrees, 

Ao  C 

It  would  be  prudent  to  evaluate  the  case  where  p  F.„  Is  stable  at 

AS 

01  .  Following  the  path  Indicated  In  Figure  5,  see  that  oi,  is  probably 
not  greater  than  oi^,  using  the  results  of  Reference  9.  In  either 
case.  It  Is,  therefore,  concluded  that  a  large  amplitude  PIO  Is  pos¬ 
sible  at  frequencies  between  oi^  -  3.3  and  oi^  «  4.18. 

Using  the  same  procedure,  the  modified  flight  control  system  (as 
defined  In  Reference  9)  would  be  PlO-free  provided  the  amplitude  of 
Fas  (when  F^g  is  approximately  sinusoidal)  is  less  than  about 
14  pounds.  It  Is  impossible  to  say  with  complete  assurance  that  this 
will  always  be  the  case.  In  comparison  with  the  prototype  system.  It 
is  nevertheless  clear  that  PIO  Is  much  less  likely  with  the  modified 
system. 

Observe  that  the  present  method  of  PIO  assessment  provides  a 
specific  frequency  estimate  for  which  closed  loop  stability  is  neces¬ 
sary  if  PIO  is  to  be  (probably)  avoided.  With  real-world  nonlineari¬ 
ties,  It  Is  very  often  not  possible  to  eliminate  the  mathematical 
possibility  that  PIO  can  occur.  The  a  priori  estimate  of  required 
closed  loop  bandwidth  provides  a  means  of  discriminating  between  the 
mathematically  possible  and  the  physically  probable  PIO  states. 
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Figure  8.  Time  History  of  an  X-15 
Controllability  Limit  (: 
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The  airplane  dynamics  are  considered  to  be  linear  in  this  flight 
condition.  The  rudder  control  was  purposely  not  used. 

First »  examine  Che  FIO  tendencies  predicted  from  closed  loop 
dynamics.  From  Figure  9,  compute 

u)  -  2.15 
c 

)t  ^  (jw^)  ■  -150  degrees 
a 

This  Indicates  no  PIO  problem  due  to  5  dynamics,  per  se. 

This  situation  changes  when  the  open  loop  dynamics  are  con¬ 
sidered.  The  dutch  roll  damping  ratio  and  frequency  ace: 

.034 
-  3.11 

Also,  >  ((>/6^(ju)j)  -  -227  degrees.  Thus,  by  the  flow  chart,  <  .2, 
and  the  possibility  for  PIO,  therefore,  exists  at  frequencies  approxi¬ 
mately  equal  to  the  dutch  roll  frequency.  In  fact,  ()>  •+  Is  unstable 
at  (0^  and  a  PIO  Is,  therefore,  predicted. 

Note  that  for  this  case  ^  The  closed  loop  Instability 

could.  In  theory,  be  arrested  with  sufficient  pilot  gain.  This  would 
require  a  gain  Increase  by  the  pilot  of  as  much  as  26  db  (a  factor  of 
19.953).  This  Is  not  a  realistic  possibility. 

C.  THE  M2-F2 

Lateral-directional  PIO  problems  with  this  aircraft  occurred 
during  three  of  16  flights  with  three  different  pilots.  These  experi¬ 
ences  are  documented  In  Reference  11.  The  aerodynamic,  mass,  and 
Inertial  data  used  for  the  present  report  are  from  Reference  12. 
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The  M2-F2  automatic  flight  control  system  consisted  of  a  simple 
SAS  with  yaw  rate  and  roll  rate  feedback  through  washout  circuits. 

The  SAS  authority  was  ±5  degrees  with  actuator  rate  limits  of 
30  degrees  per  second  In  roll.  Because  of  considerable  adverse 
aileron  yaw,  the  roll  subsidence  and  roll  spiral  modes  were  coupled  to 
produce  a  low  frequency  oscillatory  mode.  Aerodynamic  nonlinear It les 
destabilized  this  mode  as  the  trim  angle  of  attack  approached  0  degrees. 
Because  of  this,  an  aileron-rudder  Interconnect  (ARl)  was  used  to 
counter  the  aileron  yaw.  The  ARI  ratio  was  pilot-selectable  from  the 
cockpit. 

The  SAS  gains  and  the  ARI  ratio  were  different  on  each  of  the 
three  flights  for  which  PIO  was  obtained  (flights  1,  10,  and  16). 

Each  of  these  three  conditions,  plus  a  condition  for  which  PIO  was 
definitely  not  a  problem,  will  be  examined  below  as  a  separate  PIO 
case  history.  The  M2-F3  will  also  be  considered. 

In  view  of  remarks  made  elsewhere  in  this  report,  it  is  note¬ 
worthy  that  all  PIO  experiences  with  the  M2-F2  were  in  flight  phase 
Category  B.  It  is  shown  here,  however,  that  these  may  be  diagnosed 
using  criteria  derived  from  Category  A  requirements. 

M2-F2  Flight  1 

The  PIO  time  history  for  this  case  is  shown  in  Figure  10  (from 
Reference  11).  The  required  transfer  functions  were  computed,  using 
data  of  References  11  and  12,  to  be: 


A^j^(8)  »  .949  (.591)  [.498,  .171]  [.090,  4.585]  (4.229) 
^  (s)  -  .631  (.572)(.586)  [.688,  19.24]/Aj,j^(s)  (deg/deg) 


-6.280  (.511)(.587)  [-.424,  1.128]  (3.689)/Aj,j^(s) 
(deg/sec/lnch) 
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(s)  -  14.878  (.571)(.571)(-5.303)(5.60l)/A-- (s) 

CL 

(deg/deg) 

The  shorthand  notation  used  here  is: 

(l)--l 

fc  ,0)  1  ■  8^  +  2c  01  8  +  0)^ 

L  n  nj  n  n  n 

Note  that  In  the  following  three  figures,  6  Is  the  lateral  stick 

Ld 

deflection  denoted  by  elsewhere  in  this  report.  This  was  done  to 
preserve  the  notation  of  Reference  11. 

From  the  bank  angle  transfer  function,  we  can  compute 

0)  »  4.91,  ):  (jo)  >  “  -246  degrees 

AS 

Also,  a  resonant  mode  exists  at  the  open  loop  dutch  roll: 

-  .09  <  .2 

0),  »  4.585 
a 

y  (joij)  “  -211  degrees 
AS 

It  is,  therefore,  concluded  that  PIO  is  possible  with  these  dynamics. 
The  initial  PIO  frequency  will  occur  between  4.59  and  4.91  radians/ 
second.  Limit  cycle  properties  will  determine  the  resulting  steady 
state  PIO  behavior.  An  example  calculation  of  limit  cycles  is 
deferred  until  the  discussion  of  flight  16. 

The  time  history  (Figure  10)  Indicates  that  *  5.2  prior  to  PIO 
development.  This  is  in  good  agreement  with  the  predicted  value. 
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Figure  10. 


Pilot's  Input  and  vehicle  response 


Time  History  of  Pilot-Induced  Lateral-Directional 
Oscillation  on  H2-F2  Flight  1.  M  0.48;  h  >  2830  m 
(9275  ft)  to  1678  m  (5500  ft);  K  -  0.6;  K  -  0.6 
(from  Reference  11,  Figure  8)  ^ 


Throughout  the  M2-F2  FIO  occurrences,  the  pilots  remarked  about 
the  tendency  for  rudder  usage  to  promote  loss  of  control.  The 
transfer  function  suggests  that  this  is  due  in  part  to  complex,  non- 
minimum  phase  zeros  at  the  low  frequency.  Any  attempt  by  the  pilot  to 
control  dutch»roll- Induced  yaw  rate  with  rudder  would  produce  an 
immediate  loss  of  control  with  the  dynamics  shown.  The  time  histories 
suggest,  also,  that  rudder  was  used  for  control  of  p  and  (p.  For  the 
dynamics  shown,  the  linear  system  with  ^  -*■  5  will  be  highly  resonant 
at  frequencies  between  about  A.l  and  A. 6  rad Ians/ second .  If  the  pilot 
attempts  to  correct  bank  errors  at  ••  A. 91  with  5^,  he  will  Induce 
closed  loop  instability. 

The  M2-F2  Flight  10 

Figure  11  shows  the  PIO  time  history  for  this  case.  The  required 
transfer  functions  are: 

A^,j^(s)  -  .970  (.252)(.A78)  [.053,  .889]  [.856,  A. 262] 

^  (s)  -  lA.llO  (.571)(.665)  [.516,  2.356]/A^j^  (deg/deg) 

AS 

7^  (a)  -  -5.230  (.572)(.53A)  [-.626,  1.160]  (5.210)/A^, 

Arp  cl 

(deg/sec/deg) 

^  (s)  -  8.805  [1.000,  .571]  (-8.770) (8.665) /W_,  (deg/deg) 

°RP 


Compute : 


(D  -  3.18 
c 

>  (ju)  )  -  -138  degrees 

®  a  o  C 
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Pilot's  input  and  vehicle  response 


Figure  11 


.  Time  History  of  Pilot-Induced  Lateral-Directional 
Oscillation  on  M2-F2  Flight  10.  M  -  0.61;  h  ■ 
7020  m  (23,000  ft)  to  5800  m  (19,000  ft);  Kp  - 
0.4;  Kj  ■  0.6  Except  as  Noted;  Kx  “  0.49  (from 
Reference  11,  Figure  10) 


To  Illustrate  the  method,  let  us  assume  that  an  analysis  Indicates  a 
stable  limit  cycle  at  frequency  *  1*53  (the  PIO  frequency  in  the 
steady  state  condition).  Then  we  trauld  conclude  that: 

•  The  6^  loop  is  stable  in  the  linear  range  of  flight 
control  system  operation. 

•  The  predicted  handling  qualities  (based  on  the  Reference  5 
method)  are  excellent. 

•  PIO  is  not  a  problem  provided  that  the  flight  control  system 
is  not  saturated. 

e  Any  attempt  to  use  rudder  to  stabilize  dutch  roll  oscilla¬ 
tions  will  result  in  loss  of  control. 

But,  according  to  Figure  5,  since  ai  is  greater  than  u  ,  there  is  a 
possibility  that  the  limit  cycle  could  be  excited. 

This  was  done  (unintentionally)  in  flight  10  by  turning  off  the 
SAS  (see  Figure  11)  and  applying  an  aileron  pulse  at  t  >  4.25  seconds. 
The  aircraft  was  rolling  left  slightly  when  the  aileron  pulse  was 
applied.  As  a  result,  it  appeared  to  diverge.  The  pilot  corrected 
with  full  right  stick  and  quickly  turned  on  the  SAS.  It  appears  that 
after  hitting  the  stops  he  Immediately  applied  rudder  to  aid  in 
stopping  the  left  roll.  He  appeared  to  continue  to  use  rudder  plus 
aileron  to  control  roll  angle;  this  resulted  in  closed  loop  instabil¬ 
ity,  full  saturation  of  the  SAS,  and  entry  into  the  limit  cycle. 

The  aircraft  mode  corresponding  to  the  coupled  roll-spiral  is 
resonant  (damping  ratio  <  .2).  However,  at  the  modal  frequency, 
compute 

):  ■  -146  degrees 

*AS 
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The  linear  system  would  be  stable  at  this  bandwidth.  Since  the  fre¬ 
quency  Is  probably  less  than  possible  limit  cycle  frequencies,  then 
PIO  due  to  stimulation  of  the  coupled  roll-spiral  mode  by  any  means  Is 
unlikely.  (In  an  actual  design  study,  it  would  be  necessary  to  per¬ 
form  substantial  analyses  to  verify  what  Is  merely  assumed  here) . 

Therefore,  it  may  be  concluded  that,  with  a  certain  amount  of 
ingenuity,  the  possibility  of  PIO  could  have  been  forecast.  This  is  a 
particularly  interesting  example  since  the  basic  roll-to-aileron 
d)mamics  appear  to  be  very  good.  The  PIO  occurred  entirely  due  to 
flight  control  system  saturation  effects.  With  more  control  and  SAS 
authority,  there  probably  would  have  been  no  PIO  with  these  dynamics. 

M2-F2  Flight  16 

Figure  12  shows  the  PIO  time  history  for  this  case.  The  transfer 
functions  of  interest  are: 

A^j^(s)  =  .949  (.302)(.546)  [-.173,  .881]  [.667,  4.219] 

^  (s)  -  13.886  (.571)(.653)  [.491,  1.909]/A„(s)  (deg/deg) 
*AS 

~  (s)  =  -4.962  (.555)(.571)(3.999)  [-.642,  1.378]/A  fs) 
(deg/ sec/ inch) 

^  (s)  =  9.000  (.571)(.572)(-8.191)(8.084)/A_, (s)  (deg/deg) 

OrP 

Note  that  the  coupled  roll-spiral  mode  is  unstable. 
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Figure  12. 


Study  of  Figure  12  Indicates: 

•  A  (nearly)  2  Hz  oscillation  existed  in  aileron  stick  deflec¬ 
tion  in  the  pre-PIO  condition.  This  appears  to  be  pilot¬ 
generated  and  may  originate  from  the  dutch  roll  mode  with 
high  pilot  gain. 

•  The  PIO  appears  to  really  start  at  t  =  7  following  stick 
bottoming.  The  SAS  doesn't  appear  to  be  saturated  at  that 
point.  SAS  rate  saturation  occurs  at  about  7.5  with  SAS 
position  saturation  at  t  *  10.5. 

•  Rudder  pedal  is  approximately  proportional  to  bank  angle 
throughout  much  of  the  PIO. 

Compute: 

w  =  4.06,  7^  (jo)  )  =  -127  degrees 

®AS 

With  this  phase  angle,  the  pilot  rating  predicted  from  use  of  Figure  4 
indicates  that  pilot  rating  should  be  very  good  (about  3,5).  There  is 
no  known  way  to  predict  the  quantitative  effect  on  rating  due  to  the 
unstable  roll-spiral  mode.  However,  the  tlme-to-double  for  the 
unstable  mode  is  4.55  seconds.  This  probably  isn't  short  enough  to 
seriously  degrade  the  handling  qualities.  However,  any  periods  of 
pilot  inattention  to  bank  angle  control  will  result  in  roll  oscilla¬ 
tions  at  frequencies  approximately  equal  to  that  at  which  the  closed 
loop  root  locus  originating  at  this  mode  crosses  into  the  left-half 
plane.  This  frequency  is  1,43  rad Ians/ second.  The  closed  loop  root 
locus  of  the  linear  <(>  ->■  6^  system  with  a  pure  gain  pilot  is  shown  in 
Figure  13.  The  Bode  is  shown  in  Figure  14. 

The  manner  in  which  the  pilot's  limited  control  authority  could 
have  contributed  to  the  initiation  of  large  amplitude  PIO  will  be 


Figure  13.  Locus  of  Pilot-Aircraft  System  Closed  Loop 
Roots:  -*■  6^  (M2-F2  Flight  16) 

illustrated  with  the  following  example.  This  is  an  oversimplification 
of  the  actual  PIO  mechanics  since  SAS  saturation  in  position  or  rate 
isn't  considered.  The  example  does  Illustrate  how  SAS  saturation  may 
lave  been  promoted  even  though,  by  all  indications,  the  linear 

AS 

dynamics  were  very  good. 


or 


-  -  »4) 

The  sinusoidal  describing  function  for  simple  saturation  may  be 
found  In  Reference  13.  The  two  functions  are  plotted  on  the  gain 
phase  plot  of  Figure  15.  The  intersection  of  L(ja))  and  -1/N  corre¬ 
sponds  to  a  stable  limit  cycle  with  a  frequency  of  1.43  radians/ second 
and  a  control  stick  force  amplitude  of  12  pounds.  The  linear  function 
shown  is  plotted  for  a  pilot  gain  of  =  .135  pounds/degree;  for 
other  values  L(jtu)  is  shifted  vertically.  For  the  simple  saturation 
shown,  the  limit  cycle  frequency  is  constant.  For  a  limit  cycle  to 
exist  at  all,  an  approximately  periodic  input  control  force  of  ampli¬ 
tude  greater  than  about  10  pounds  is  required.  Since  the  limit  cycle 
is  stable,  saturated  control  oscillations  at  frequencies  greater  than 
(jjj^  =  1 .43  will  eventually  decay  to  the  steady  state  limit  cycle 
condition. 


•1/N 


Figure  15.  Gain  Phase  Plot:  M2-F2  Flight  16  (Simple  Saturation) 
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Now,  consider  the  implications  of  these  results  to  the  assessment 
of  PIO.  Following  the  procedure  of  Figure  5: 

1.  The  linear  system's  crossover  frequency  to^  -  4.06;  this  Is 
greater  than  the  limit  cycle  frequency  u  >  1.43. 

Li 

2.  The  unstable  roll-spiral  mode  almost  guarantees  that  closed 
loop  oscillations  will  be  experienced  In  flight.  Thus, 
periodic  F  is  inevitable. 

3.  The  stick  force  required  for  saturation  of  stick  deflection 
is  only  about  10  pounds.  Full  roll  control  yields  only 
about  7  degrees/ second  of  roll  rate. 

4.  The  sluggish  roll  response,  combined  with  the  likelihood 
that  closed  loop  oscillations  will  occur  originating  from 
the  roll-spiral  mode  suggests  that  the  conditions  required 
for  the  stable  limit  cycle  will  be  obtained. 

Consider  the  case  where  the  limit  cycle  conditions  are  marginally 

met;  w  =  1.43  and  F  amplitude  -  10  pounds. 

L  Ad 

1-^1  -  6.4  db  -  2.09  deg/deg 
°AS 

The  amplitude  of  6  ,  referenced  to  6  ,  Is  10  degrees  (which  corre- 

Am  3 

sponds  to  F.-  ■  10  pounds).  Then  ^  ■  20.9  degrees  and  P  - 

30  degrees/ second.  The  maximum  Increment  of  aileron  deflection  due  to 
the  roll  rate  feedback  (gain  ■■  0.2  degree/degree)  Is  6.0  degrees. 

Since  the  SAS  authority  limit  Is  only  5  degrees,  then  It  appears  that 
saturation  of  6  will  also  saturate  the  SAS.  Thus,  the  limit  cycle 

Am 

condition  discussed  above  Is  an  oversimplification;  It  Illustrates  the 
approach,  however,  which  was  the  Intent.  In  design  practice,  this 
would  be  sufficient  Information  to  warrant  a  design  review.  A 
possible  remedy  Is  to  Increase  the  pilot  control  authority;  this.  In 
fact,  was  done  for  the  M2-F3. 


t 
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The  use  of  rudder  to  augment  aileron  roll  authority  vrould  result 
In  closed  loop  Instability  and  SAS  saturation. 

It  is  believed  that  the  above  analyses  are  in  reasonable  agree¬ 
ment  with  the  actual  PIO  experience  on  flight  16. 

M2-F2  (Non-PIO  Condition) 


For  this  example,  the  aircraft  configuration  is  identical  to  that 
of  flight  16  except  for  the  trim  angle  of  attack.  The  value  used  was 
oq  =  8.1  degrees  (it  was  -2.3  degrees  for  the  flight  16  PIO  condition). 

For  this  case: 


Acl(s)  =  .949  (.581)(1.152)  [.851,  .222]  [.202,  8.572] 

^  (s)  -  7.129  (.571)(.583)  [.207,  4.435]Mpj^(s)  (deg/deg) 


w  =  3.31 
c 

)  7^  (jw  )  =  -138  degrees 
A  c 


Observe  that  the  equivalent  <  1.  Because  of  this,  there  are 

no  resonant  closed  loop  modes.  The  damping  ratio  for  the  dutch  roll 
mode  is  close  to  the  0.2  criterion  for  open  loop  resonance.  However, 
u)d  >  by  a  substantial  margin.  It  is  concluded  that  neither  closed 
nor  open  loop  initiated  PIO  are  likely  for  this  configuration.  Note 
that,  even  with  saturated  6^,  there  is  no  resonance  to  excite  PIO. 


The  rudder  dynamics  for  this  case  (not  shown)  permit  normal 
rudder  usage,  in  contrast  to  the  situation  at  low  angle  of  attack. 


D.  M2-F3  (SAS  ON) 


The  M2-F2  was  modified  by  adding  a  fixed  center-fin.  This 
decreased  the  adverse  aileron  yaw  and  eliminated  the  requirement  for 
the  ARI.  Control  authority  was  doubled.  The  example,  here,  is  for 
the  same  flight  condition  as  the  PIO  case  for  flight  16  of  the  M2-F2. 

Aj,j^(s)  -  .949  [.157,  .853]  (.992,  .456]  [.685,  3.876] 

^  (s)  -  16.496  (.571)(.612)  [.343,  2.8581/A^j^(s)  (deg/deg) 

(s)  -  -4.962  (.547)(.571)  [-.629,  1.339]  (4. 296) /A.,  (s) 
(deg/sec/deg) 

^  (s)  =•  9.000  (-8.191)(.571)(.572)(8.084)  (deg/deg) 

HP 

Li  -  3.05 
c 

y.  (joj  )  «  -150  degrees 
AS 

As  a  first  approximation,  it  is  reasonable  to  expect  that  limit  cycles 
will  exist  at  frequencies  equal  to  those  for  the  original  M2-F2.  The 
Increased  control  authority  will  reduce  the  probability  that  motion 
amplitudes  sufficient  to  produce  a  limit  cycle  will  be  experienced  in 
flight. 

Exercise  of  the  criteria  from  Figure  5  suggests  that  PIO  is 
unlikely  due  to  aileron  control.  However,  use  of  rudder  for  control 
of  bank  or  yaw  rate  can  destabilize  the  system  with  closed  loop  band¬ 
width  about  equal  to  1  rad Ian/ second.  This  could  lead  to  control 
system  saturation  and  fully  developed  PIO.  A  complete,  detailed 
analysis  of  the  nonlinear  system  would  be  required  to  establish  the 
conditions  by  which  PIO  could  be  sustained  due  to  rudder  usage. 
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E.  M2-F3  (SAS  OFF) 


!l 


This  case  was  simulated  but, 
trim  ag  •  -4  degrees. 

A^^(s)  »  .949  [-.575,  . 

7^  (s)  =  16.496  [.159, 
^AS 

~  (s)  -  -4.962  [-.986 
^  (deg/sec/deg) 


apparently,  not  flight  tested. 

812]  [.475,  3.251] 

2.958]/A^j^(s)  (deg/deg) 

,  1.425]  (3.629)/Aj,j^(8) 


The 


^  (s)  =  9.000  (-8. 191) (8. 084)  (deg/deg) 
'^RP 


w  -  3.39 
c 

y  (jw  )  =  -164  degrees 
AS 


These  dynamics  are  similar  to  those  for  the  M2-F2  flight  16  PIO 
case.  A  detailed  analysis  would  probably  Indicate  that  limit  cycles 
are  unlikely  due  to  the  increased  control  authority. 

The  rudder  dynamics  are  poor.  Closed  loop  control  of  ([i  or  r  with 

6  will  produce  closed  loop  instability.  However,  control  can  prob- 
Kr 

ably  be  recovered  by: 

•  ceasing  rudder  usage,  and 

•  closing  the  ^  6^  loop  in  a  normal  manner. 

These  analytical  observations  are  consistent  with  the  simulator 
tests  of  this  configuration.  Figure  16  shows  a  system  instability 
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t.  aac 


Figure  16.  M2-F3  S lx -Degre.e-of -Freedom  Simulator  Time 

History  of  Final  Turn  and  Approach  to  Landing 
with  SAS  Off  (from  Reference  11,  Figure  23) 


induced  by  attempted  rudder  coordination  with  aileron.  A  roll  diver¬ 
gence  was  created;  recovery  was  accomplished  in  one  cycle,  without 
control  saturation.  The  frequency  of  the  large  amplitude  oscillation 
Is  approximately  equal  to  the  frequency  at  which  the  ^  6  root 

d 

locus,  originating  from  the  unstable  roll-spiral  mode,  crosses  into 
the  left-half  plane. 


F.  ESP  CONFIGURATIONS 

The  Equivalent  Systems  Program  (ESP)  data  was  introduced  in 
Section  II.  The  data  source  was  Reference  8.  The  analysis  of  these 
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1 

i 

data  according  to  the  requirements  of  Figure  5  was  particularly  easy. 

The  simulated  aircraft  dynamics  are  all  linear.  Only  the  aileron 
control  was  used.  Further,  the  dutch  roll  mode  was  canceled  from  the 
roll  transfer  function.  The  roll  dynamics  consisted  of  conventional 
spiral  and  roll  subsidence  modes  plus,  for  some  configurations,  a  low 
pass  control  stick  output  filter  and  substantial  time  delay. 

The  computed  data  are  shown  in  Table  1.  Since  there  are  neither 
limit  cycles  nor  resonant  open  loop  dynamics  to  consider,  only  the 
crossover  conditions  are  required  for  testing  the  proposed  PIO  analy¬ 
sis  method.  The  PIO  frequencies  are  shown  for  the  four  cases  analyzed 
at  McDonnell  Aircraft  (Reference  8);  no  time  history  data  were  avail¬ 
able  for  the  other  cases.  The  Cooper -Harper  ratings  are  shown. 

Including  replications. 

From  Table  1  It  may  be  concluded  that  the  simple  phase  criterion, 
evaluated  at  yields  very  encouraging  results.  Only  configura¬ 
tion  L7A  shows  a  clear  contradiction  between  theory  and  test, PIO 
predicted  but  not  experienced.  However,  L7A  was  only  flown  once. 

For  several  other  configurations  for  which  PIO  was  predicted  to 
be  a  problem,  but  was  not  found  in  test,  the  pilot  comments  were 
consistent  with  those  features  that  could  produce  PIO  in  less  favor¬ 
able  circumstances.  Representative  comments  are  shown  in  Table  1. 

Even  if  these  cases  prove  eventually  to  be  not  susceptible  to  PIO,  the 
phase  criterion  is  still  useful  since  it  does  indicate  configurations 
which  have  possible  handling  qualities  deficiencies,  even  though  not 
of  the  PIO  variety. 
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TABLE  1.  ESP  FLIGHT  TEST:  PIO  ASSESSMENT 


Conf 

00 

c 

AS 

Predic¬ 
tion:  PIO 
Possible? 

Flight 

Test: 

PIO 

Obtained 

“pio 

Pilot 

Rating 

LI 

3.99 

-168 

no 

no 

4 

L2 

3.70 

-179 

marginal 

no 

3 

L3 

4.04 

-164 

no 

no 

4 

L4 

3.63 

-176 

no 

no 

4 

L4A 

3.63 

-176 

no 

no 

3 

L5 

3.53 

-169 

no 

no 

2/2/2 

L6 

3.67 

-170 

no 

no 

2 

L7 

3.44 

-183 

yes 

no 

3 

L7A 

3.44 

-201 

yes 

no 

4 

L8 

2.67 

-204 

yes 

yes 

2.3-3.14 

5 

L8A 

2.58 

-208 

yes 

(1) 

6 

L8B 

2.48 

-212 

yes 

yes 

9/5 

L9 

3,77 

-181 

marg Inal 

no 

2 

LlOA 

3.77 

-192 

yes 

yes 

5 

Lll/A/B 

3.77 

-205 

yes 

(2) 

3/4/5/6 

LllC 

3.77 

-226 

yes 

yes 

4.8 

9/6 

LllD 

2.77 

-225 

yea 

yes 

10 

L12 

3.25 

-179 

marginal 

no 

4/5 

L13 

3.16 

-186 

yes 

(3) 

4 

L14 

3.02 

-194 

yes 

(4) 

5/7 

LI4A 

2.48 

-210 

yes 

yes 

8 

L14B 

2.12 

-222 

yes 

yes 

2.09 

8/10 

L15 

3.34 

-192 

yes 

(5) 

4/5 

LI  6 

3.34 

-202 

yes 

(6) 

3/4 

L16A 

3.34 

-214 

yes 

yes 

3.5 

8 

Comments: 

I 

I 

I  (1)  Overshoots  bank  angle  with  3-A  slow  oscillations. 

‘  (2)  Rated  sluggish  and  sensitive  in  bank  angle  control;  some  require¬ 

ment  for  opposite  6  to  stop  roll  rate;  slight  overcontrol 
tendenc les. 

(3)  Oscillated  about  final  bank  angle;  tended  to  overturn  in  sidesteps. 

(4)  Overshoot  desired  bank  angle,  then  oscillate;  could  not  bring 
myself  to  put  in  large  corrections  during  sidestep. 

(5)  Sidestep  forced  bang-bang  control;  overcontrolled  bank  angle. 

(6)  Overcontrol  in  sidestep. 
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G.  REENTRY  VEHICLES  (HARPER) 


Simulation  of  reentry  vehicle  lateral-directional  dynamics  in  a 
variable  stability  aircraft  was  reported  in  Reference  14.  One  hundred 
and  twenty  nine  configurations  were  tested.  The  piloting  tasks  were 
of  the  flight  phase  Category  B  type,  plus  IFR  tracking.  One  evalua¬ 
tion  pilot  was  used;  only  three  repeat  runs  were  made.  The  primary 
purpose  of  the  test  was  to  accumulate  data  for  the  general  evaluation 
of  lateral-directional  handling  quality  parameters  such  as  and 

PIO  was  not  specifically  sought  and,  in  fact,  was  only  found 
in  two  cases  (configurations  94  and  102).  For  configuration  90, 
undanped  roll  oscillations  were  found  during  IFR  tracking  which  were 
not  a  problem  in  VFR  control.  Consider  each  of  these  three  configura-- 
tions: 


1.  Configuration  90: 


0)  =  3.58  >  (jw  )  =  -160 

^AS 


No  PIO  problem  is  indicated  in  closed  loop  roll  control.  The  dutch 
roll  mode  damping  ratio,  is  .268;  the  mode  is  not  resonant 
according  to  the  definition  used  here.  Note,  however,  that 

=  -179.  It  is  not  inconceivable  that  due  to  instrument  lags  or 
limited  (^-resolution,  either  the  open  loop  dutch  roll  oscillations  did 
create  an  IFR  control  problem  or  the  resultant  phase  angle  at 
violated  the  -180  degrees  criterion.  This  is  conjectural,  but  not  an 
unreasonable  account  of  the  differences  observed  between  VFR  and  IFR 
flight. 

2.  Configuration  94: 

to  =  3.43  ):  (ju)  )  »  -149 

C  0  c 

°AS 
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No  PIO  problem  is  indicated  at  the  predicted  crossover  frequency.  The 
dutch  roll  mode  is  resonant  ■  .04). 

0),  =  4.55  )  7^  *  -230  degrees 

^AS  ** 

This  indicates  a  possible  PIO  problem  at  the  dutch  roll  frequency. 

3.  Configuration  102: 

The  results  of  this  case  are  very  similar  to  the  results  for 
configuration  94: 

u  =  2.88  )  7^ —  (ju)  )  =  -147 

^AS 

^d  =  *15 

w  =  4.15  ^  (jo).)  *  -215 

AS 

According  to  the  proposed  theory,  there  are  several  other  con¬ 
figurations  that  would  have  been  designated  as  possibly  PlO-prone. 
These  are:  54,  77,  78,  85,  86,  92,  93,  101,  and  105.  These  were  not 
identified  as  PlO-prone  in  flight  test.  Whether  this  was  due  to 
insufficient  testing,  to  test  technique,  or  to  a  genuine  lack  of  PIO 
problems  cannot  be  determined.  It  is  interesting  that  for  all  these  9 
configurations,  except  92,  >  1.  The  pilot  ratings  were  6-9  for 

all  cases  except  54,  85,  and  105. 

It  is  tentatively  concluded  that  the  PIO  analysis  method  of 
Figure  5  gives  realistic  Indications  of  PIO  problems,  specifically, 
and  potential  handling  problems  for  most  of  the  129  cases  tested. 
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H.  THE  PRINCETON  DATA 


Reference  15  documents  flight  tests  of  134  lateral-directional 
dynamic  configurations  simulated  with  the  Princeton  Navion.  Perhaps 
what  this  report  shows  more  than  anything  else  is  the  total  breakdown 
of  the  empirical  approach  to  handling  qualities  when  applied  to  the 
study  of  lateral-directional  dynamics.  There  are  too  many  parameters 
and  their  interactions  are  far  too  complex  for  any  practical  test 
program  to  deal  with  in  a  strictly  empirical  manner.  The  authors 
noted,  for  example,  the  inability  to  independently  vary  and  |(j>/3| 
over  a  useful  dynamic  range.  One  Interesting  case  for  PIO  research 
was  encountered. 

Configuration  81  was  found  to  have  fair  handling  qualities  in  the 
approach  to  carrier  landing  task  (R  =  4.5,  Cooper  scale).  No  actual 
touchdowns  were  made,  nor  was  flare  control  a  requirement.  However, 
this  same  configuration  demonstrated  PIO  problems  in  nontracking 
flight  tasks  (Category  B) .  For  the  configuration  81  dynamics: 

u  =3.24 
c 

^  6^  (j“c)  “  -^51 

AS 


^d  = 


u)j  =  1*8 

a 


>  (J«j)  -  -199 

AS 


These  data  are  consistent  with  the  flight  test  results.  No  PIO  prob¬ 
lems  in  tracking  are  Indicated.  The  dutch  roll  mode  is  resonant,  and 
a  possible  PIO  problem  at  is  indicated. 


* 
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I. 


CONCLUSIONS 


Almost  all  of  the  data  considered  support  the  proposed  PIO 
assessment  method  of  Section  IV.  In  those  cases  where  the  method 
Indicated  possible  FIO  problems  to  occur,  but  where  none  were  found  In 
flight  test,  the  handling  qualities  were  nearly  always  poor, 
nonetheless. 


SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 


Single  axis,  roll-only  PIO  appears  to  be  possible.  That  is, 
lateral-directional  PIO  does  not  necessarily  require  simultaneous 
pilot  control  of  multiple  cues  (e.g.,  4i  and  r)  or  the  use  of  more  than 
one  control.  While  PIO  may,  in  some  cases,  be  initiated  through  such 
mechanisms,  the  data  used  in  this  study  indicate  that  possible  PIO 

problems  can  generally  be  forecast  based  on  (|)  ->  F._  dynamics. 

Ao 

The  method  proposed  for  the  assessment  of  lateral-directional  PIO 
is  very  similar  to  traditional  methods.  It  differs  mainly  in  the 
quantitative  prescription  of  specific  frequencies  for  exercise  of 
stability  criteria  (linear  or  limit  cycle). 

It  isn't  clear  whether  lateral-directional  PIO  is  inherently 
"different"  from  longitudinal  PIO  or  not.  The  longitudinal  theory 
(Reference  1)  hypothesized  a  coupling  between  the  pilot's  control  of 
pitch  attitude  and  normal  acceleration  as  the  source  for  PIO  diffi¬ 
culties.  During  the  present  work,  an  attempt  was  made  to  extend  the 
longitudinal  theory  to  the  lateral-directional  case.  There  is  some 
weak  evidence  that  lateral  acceleration  may  be  a  factor  in  some 
lateral-directional  PIO.  Unfortunately,  lateral  acceleration  data  are 
almost  never  reported  in  the  handling  qualities  literature.  Because 
of  this,  it  wasn't  possible  to  pursue  this  explanation  for  PIO. 
However,  in  the  spirit  that  the  simplest  theory  is  the  best,  the 
roll-only  PIO  model  may  be  satisfactory  until  new  data  suggests 
otherwise. 

One  area  where  lateral  acceleration  may  play  a  role  is  with 
direct  side  force  control  applications.  There  are  no  data  known  to 
this  author  about  PIO  with  such  flight  control  systems.  Still,  it  is 
reasonable  to  assume  that  PIO  could  occur  in  tasks  such  as  ground 
attack  or  level  turns  for  bombing.  If  so,  the  longitudinal  theory  may 
®PPly.  directly,  provided  lateral  acceleration  and  an  appropriate 


60 


outer  loop  (guidance)  cue—whlch  is  task  dependent — are  used  In  place 
of  normal  acceleration  and  pitch  attitude. 

The  apparent  success  of  the  roll^only  PIO  model  suggests  that  the 
pitch  PIO  theory  should  be  revised  to  account  for  a  single  axis  pitch 
attitude  PIO.  The  assessment  method  would  be  Identical  with  that 
developed  In  this  report  for  roll  PIO.  A  corollary  to  this  single 
axis  pitch  PIO  mode  Is  that  PIO  can  be  obtained  In  fixed  base  simula¬ 
tion  which  Is  physically  equivalent  to  that  found  In  flight  test.  It 
Is  noted,  however,  that  In  the  absence  of  system  nonlinear It les, 
flight  control  system  time  delays,  or  higher  order  system  dynamics, 
the  stability  criterion  ):  0/F„„(jw  )  >  -180  degrees  will  almost 
always  be  satisfied;  thus  single  axis  pitch  PIO  probably  can't  occur 
unless  one  or  more  of  these  effects  are  present.  This  Is  not  the  case 
with  the  original  theory  of  Reference  1.  In  other  words,  the  classic 
airplane  can,  according  to  Reference  1,  exhibit  PIO  due  to  coupling, 
through  pilot  control,  of  pitch  and  normal  acceleration.  However,  by 
the  theory  of  this  report.  It  would  be  unlikely  to  have  a  pitch-only 
PIO. 

It  Is  recommended  that  future  flight  tests,  conducted  for  the 
purpose  of  extending  the  PIO  data  base  (longitudinal  or  lateral- 
directional)  ,  should  be  done  using  direct  force  control  techniques  to 
decouple  attitude  and  linear  acceleration  modes.  It  should  be  pos¬ 
sible  to  Identify  the  various  effects  on  PIO  or  handling  qualities 
when  the  testing  Is  done  In  this  way.  There  Is  a  certain  urgency  to 
accumulate  data  of  this  sort  In  view  of  the  flight  control  system 
state  of  the  art. 

No  revisions  to  MIL-F-878SC  are  recommended  based  upon  analyses 
In  this  report.  More  handling  qualities  data  are  required  from  flight 
tests  dedicated  to  testing  proposed  criteria  before  such  proposals  can 
be  made — If,  Indeed,  they  are  warranted. 


APPENDIX 

HANDLING  QUALITIES  PREDICTION  IN  APPROACH  AND  LANDING 


Reference  16  documents  one.  of  the  more  recent  and  Important 
handling  qualities  data  bases  to  become  available.  Longitudinal 
dynamics  were  varied  over  a  broad  range  and  simulated  with  the  vari¬ 
able  stability  T-33.  The  piloting  task  was  landing  approach,  with 
actual  touchdowns  In  most  cases.  This  data  set  Is  Identified  by  the 
acronym  LAHOS  (Landing  Approach  Higher  Order  Systems) . 

Longitudinal  PIO  In  the  approach  and  landing  task  were  specifi¬ 
cally  sought  as  a  subset  of  the  LAHOS  experiment.  Several  were  found. 
This  data  set  Is  important  to  the  present  work  because  It  Increases 
the  general  knowledge  of  the  PIO  problem,  and  permits  an  Independent 
check  on  the  handling  quality  parameters  proposed  In  Reference  5.  Of 
the  latter,  the  parameters  0/F„{jt»  )  and  w  are  the  most  Important. 

£iw  C  C 

The  other  Is  the  normal  acceleration  phase  parameter: 

3np 

<t>  (ju)  )  a  ):  —  (jw  )  -  14.3  w 
c  F_-  c  c 

ES 

The  parameters  u  and  ^  0/F„„(jta)  )  are  of  central  Importance  to 

C  wO  C 

the  method  for  PIO  assessment  proposed  In  the  main  body  of  this 
report. 

It  Is  noted  that  these  two  parameters  have  already  been  evaluated 

for  the  LAHOS  data  set  (In  Reference  17).  The  calculation  of  u  was 

c 

made  using  the  original  formula  (Reference  5).  Here,  the  slightly 
modified  version  of  the  crossover  frequency  formula  from  Section  II  Is 
used: 


bi  -  6.0  +  0.24  X  S 
c 


uaoMoim  BUM  ujml-mut  numi 
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where 


S  «  average  slope  (in  db/octave)  of  |6/Fgg(ju>)l  on  1  s  u)  £ 

6  rad Ians/ second 

In  fact,  this  formula  was  developed,  in  part,  from  inspection  of  the 
LAHOS  data.  However,  no  attempt  was  made  to  optimize  the  u)^  formula 
in  any  sense.  The  use  of  the  modified  formula  yields  significantly 
better  correlation  between  pilot  rating  and  6/F„(.(ju  )  and  generally 
appears  to  be  consistent  with  observable  closed  loop  oscillation 
frequencies. 

This  correlation  is  shown  for  all  LAHOS  configurations  in 
Figure  A1 .  The  empirical  rating  boundaries  from  Reference  5  are 
superimposed  on  the  data  points  (these  are  the  same  boundaries  shown 
in  Figure  A,  Section  IV). 


Figure  Al.  Correlation  of  LAHOS  Data  with  ):  e/F„„(j(A)  ) 

bs  c 


with  the  exception  of  configurations  1-A,  2-A,  A-0,  and  5-1,  all 
the  other  data  points  show  consistent  variation  with  )f-  • 

The  spread  in  the  data  at  fixed  about  2  1/2  Cooper -Harper 

units — over  twice  that  of  the  same  fit  achieved  for  the  Neal-Smlth 
data  (Reference  18)  analyzed  in  Reference  5  and  shown  in  Figure  A. 

Reference  5  proposed  a  handling  qualities  criterion  based  on  the 
normal  acceleration  phase  angle  parameter  4>  as  follows: 

Level  1:  4>  2  -160  degrees  for  any  value  of  ):  6/F_„(ju)  )  S 

iio  C 

-130  degrees,  or  any  4  when  )  e/F__(Ja)  )  2  -122  degrees 

Level  2:  4  >  -220  degrees  for  any  value  of  )[  0/F  „(ju)  )  > 

£*0  C 

-165  degrees, or  any  4>  when  ):  e/Fgg(ju^)  2  -1A8  degrees 

Level  3:  otherwise 

When  these  criteria  are  applied,  it  can  be  shown  that: 

•  Configuration  2-A  is  predicted  to  be  a  definite  Level  2 
condition:  ^  6/F_„(j(i)  )  =  -133  degrees,  4>  (jo)  )  » 

c«o  C  C 

-20A  degrees. 

•  Configurations  1-A,  A-0,  and  5-1  are  all  predicted  to  be 
borderline  Level  2/3  since,  in  each  case,  ) 

within  2  degrees  of  the  borderline  value  (-1A8  degrees) . 

It,  therefore,  appears  that  these  four  cases  are  entirely  con¬ 
sistent  with  the  handling  quality  assessment  rules  developed  in  Refer¬ 
ence  5  based  on  air  combat  type  control  tasks. 

Of  those  LAHOS  configurations  which  obey  the  rating  variation 
from  Reference  5,  only  A-C  would  be  mlsrated  by  imposing  the  4>  (Ju^) 
criteria.  For  A-C,  ♦  (ju  )  ■  -225  degrees  and  >  0/F_-(ju  )  “ 

-137  degrees.  According  to  the  cited  criteria,  this  configuration 
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should  also  have  been  Level  2  (l.e. ,  similar  to  the  other  4).  This 
configuration  was  flown  twice.  In  both  Instances,  the  pilot  com¬ 
plained  of  pitch  bobble  and  some  abruptness  In  Initial  pitch  response. 

There  are  23  configurations  that  lie  significantly  below  the 
Category  A  boundaries  from  Reference  5.  There  are  several  possible 
explanations  for  this: 

1.  The  >  parameter  is  Invalid. 

2.  The  ):  parameter  does  not  apply  to  the  approach 

and  landing  task. 

3.  There  are  additional  (but  unknown)  handling  quality  metrics 
which  will  unify  the  two  data  sets. 

4.  There  were  systematic  differences  between  the  processes  by 
which  the  two  data  sets  were  collected,  and  this  resulted  in 
skewing  of  the  one  set  relative  to  the  other. 

5.  The  process  by  which  one,  or  both,  data  sets  were  collected 
was  faulty. 

Considering  the  self-consistency  of  the  LAHOS  data  when  corre¬ 
lated  with  y  6/F„„(ju)  )  [with  necessary  allowance  for  <!>  (jw  )]  and  its 
CiO  c  c 

successes  elsewhere,  the  first  two  possibilities  are  difficult  to 
accept.  The  third  possibility  Is  impossible  to  evaluate;  there  are 
undoubtedly  other  metrics  which  could  be  correlated  with  some  or  all 
of  these  data.  Time  response  calculations  with  the  LAHOS  configura¬ 
tions  suggest  that  a  rise-time  criterion  on  pitch  rate  has  some  merit, 
for  example.  The  fifth  possibility  will  be  discounted  out-of-hand. 

There  were  at  least  two  possible  effects  which  may  have  produced 
systematic  differences  between  the  Neal-Smlth  data  set  and  the  LAHOS 
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data.  The  most  apparent  difference  is  in  the  nature  of  the  tasks 
performed  by  the  evaluation  pilots.  One  was  terminal,  the  other  was 
not.  A  second,  but  related,  effect  may  have  been  approach  and  landing 
performance,  and  its  effect  on  the  handling  qualities  evaluations. 

This  is  one  of  the  grey  areas  of  present-day  handling  qualities 
research.  Qualitative  evidence  suggests  that  good  performance  in  a 
manual  control  task  can  result  in  substantially  improved  pilot 
ratings  when  compared  with  nominal  rating  data. 

It  is  not  entirely  clear  how  an  evaluation  pilot  determines 
control  performance  in  a  landing  task.  The  really  poor  configurations 
in  Category  A  tasks  probably  tend  to  be  poor  in  the  landing  task  as 
well.  This  may  account  for  the  consistency  observed  between  the  LAHOS 
data  and  the  Reference  5  boundaries  in  the  Level  3  region.  In 
Category  A  (tracking)  tasks,  pitch  attitude  provides  a  continuous 
display  of  control  performance;  there  is  no  ambiguity  about  perform¬ 
ance  of  the  control  task.  Although  there  are  physical  cause  and 
effect  relations  between  attitude  variations  and  touchdown  conditions 
in  the  power  approach,  the  relations  are  not  necessarily  clear  to  an 
evaluation  pilot,  moment-to-moment. 

The  use  of  Cooper-Harper  ratings  in  the  terminal  control  task  may 
also  be  questioned.  The  scale  is  defined  in  terms  of  "workload"  and 
"performance."  In  an  approach,  the  pilot  workload  may  be  very  low, 
even  with  degraded  handling  qualities,  provided  the  aircraft  is 
trlauned  at  an  acceptable,  nominal  state  of  motion.  In  flare,  the 
workload  may  be  intense  for  a  few  moments.  It  isn’t  clear  how  a  pilot 
can  construct  a  rational  rating  when: 

•  the  approach  control  was  nearly  nominal,  without  significant 
disturbances, 

s  a  brief  period  of  high  workload  was  experienced  in  flare  and 
postflare  due  to  unanticipated  pitch  overshoots  which 
excited  undesirable  heaving  motions. 
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with  a  few  well-timed  control  Inputs,  the  undesirable 
motions  were  arrested,  and 


•  the  aircraft  was  landed  safely,  within  nomlnals. 

Based  mainly  on  pilot  comments  from  Reference  16,  It  appears  that 
scenarios  such  as  this  were  not  uncommon  during  the  LAHOS  trials.  The 
resulting  Cooper-Harper  rating  may  very  well  depend  upon  how  much  the 
pilot  blames  himself  for  the  unanticipated  response,  or  upon  how 
smooth  his  control  inputs  were.  Because  of  considerations  of  this 
type,  it  is  conceivable  that  aircraft  having  very  poor  handling  In 
Category  A  tasks  could  be  rated  very  favorably  in  the  approach  and 
landing  task  with  the  same  dynamics.  In  such  cases,  however,  one  may 
expect  that  eventually  the  situation  will  occur  when  the  pilot  will 
excite  the  undesirable  responses. 

This  jaundiced  view  is  supported  by  extensive  fixed-base  simu¬ 
lator  results.  Aircraft  dynamics  similar  to  those  of  a  current  ground 
attack  aircraft  were  simulated.  A  terrain  board-TV  projection  visual 
display  was  used.  The  aircraft  dynamics  were  modified  by  adding  time 
delay  increments  to  the  solutions  of  the  aircraft  equations  of  motion. 
Since  time  delay  has  no  effect  on  the  computed  crossover  frequency, 
u)^,  all  of  the  pilot  rating  data  obtained  as  a  function  of  time  delay 
were  for  constant  ui^.  The  resulting  correlations  between  ratings  and 
):  ®re  shown  in  Figure  A2  for  two  tasks:  approach  and 

landing  (the  circles)  and  ground  attack  (the  filled  triangles).  All 
of  the  data  for  4  test  pilots  are  shown  for  six  values  of  time  delay. 
The  rating  boundaries  derived  in  Reference  5  are  shown  on  the  figure, 
also. 


Comparison  of  Figures  A1  and  A2  suggests  that  the  simulator  data 
for  the  approach  and  landing  task  is  qualitatively  similar  to  the 
LAHOS  data — especially  for  >  ^  ”200  degrees.  The  simulated 

ground  attack  task,  however,  was  substantially  more  difficult  than  the 
landing  task,  with  the  exception  that  for  very  large  lags  in  Y  d/F^^, 
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FIXED  BASE  SIMULATOR  DATA : 
O  APPROACH  a  LANDMQ  TASK 


Figure  A2.  Correlation  of  Fixed  Base  Simulator  Data 
with  >  0/Fgg(ju)^) 


both  tasks  were  often  rated  poorly.  Figure  A3  is  a  plot  of  the 
ratings  (averaged  at  each  value  of  delay)  versus  delay  for  the  two 
tasks  in  Figure  A2  plus  a  terrain  following  task. 


On  the  basis  of  average  pilot  ratings,  the  ground  attack  task  is 
clearly  the  most  difficult.  However,  the  ratings  appear  to  originate 
from  a  statistical  process.  On  occasion,  it  is  expected  that  ratings 
in  the  approach  and  landing  task  will  equal  the  average  rating  found 
in  the  Category  A  task — even  though  there  are  significant  differences 
between  the  average  ratings  for  the  two  tasks  when  sufficient  rating 
data  have  been  collected.  It  was  also  determined  (qualitatively)  that 
in  the  approach  and  landing  task  those  flights  for  which  the  pilot 
ratings  were  very  optimistic  were  also  those  for  which  task  perform¬ 
ance  was  very  good  with  minimum  control  activity  by  the  pilot.  For 
these  cases,  it  appears  that  the  pilot  was  able  to  establish  a 
trimmed,  nominal  approach  which  required  minimum  workload.  This  even 
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COOPER-HARPER  RATING 


TIME  DELAY*- ms 

Figure  A3.  Raw  Data:  Longitudinal  Task  Comparisons 
Averaged  Pilot  Data 
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occurred  with  configurations  that  were  basically  unflyable  In  the 
ground  attack  task.  It  would  be  Interesting  to  determine  for  the 
LAHOS  data  whether  comparable  situations  existed.  This  might  be  done 
by  spectral  analysis  of  the  pilot's  control  response  to  determine 
whether  periods  of  nontracking  control  existed. 

One  conclusion  that  may  be  drawn  from  this  comparison  Is  that  It 
may  be  unconservative  to  base  handling  quality  design  specifications 
for  approach  and  landing  on  data  obtained  from  approach  and  landing 
testing — at  least  not  without  exercising  some  very  critical  judgments. 
In  view  of  the  statistical  nature  of  the  landing  accident,  we  ought  to 
be  more  skeptical  about  the  assumption  that  "looser,"  more  liberal 
design  criteria  are  satisfactory  for  Category  C  tasks  than  would  be 
the  case  for  Category  k. 

It  Is  recommended  that  these  speculations  be  tested,  first  In  a 
moving-base  simulator  and.  If  warranted,  later  in  flight.  It  would  be 
of  great  value  to  determine  whether  the  LAHOS  configurations  would  be 
given  ratings  consistent  with  those  of  the  Neal-Smith  configurations 
when  comparable  tasks  were  flown.  It  is  also  desirable  that  the 
ability  of  the  parameters  u)  ,  >  0/F„_(ju)  )  and  ^  (ju  )  to  unify  a 
spectrum  of  precision  tracking  data  (longitudinal)  be  fully  tested. 

Until  It  Is  proven  otherwise  with  suitable  testing.  It  will 
remain  the  opinion  of  this  author  that  the  approach  and  landing  data 
of  Reference  16  generally  follows  the  variations  derived  in  Refer¬ 
ence  5  based  on  Category  A  data.  Those  data  from  Reference  16  that 
differ  substantially  from  the  boundaries  derived  In  Reference  5,  with 
y  0/Fp„(j(»)  )  as  the  principal  correlator,  reflect  the  peculiarities  of 
obtaining  subjective  opinion  ratings  In  terminal  control  tasks. 
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